Computational Design of Axion Insulators Based on 5d Spinels Compounds 
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Based on density functional calculation with LDA+U method, we propose that hypothetical Os- 
mium compounds such as CaOs2 04 and SrOs2 04 can be stabilized in the geometrically frustrated 
spinel crystal structure. They also show some exotic electronic and magnetic properties in a rea- 
sonable range of on-site Coulomb correlation U such as ferromagnetism and orbital magnetoelectric 
effect characteristic to Axion electrodynamics. Other electronic phases including 3D Dirac metal 
and Mott insulator exist and would make perspective 5d spinels ideal for applications. 



PACS numbers: 73.43.-f, 72.25.Hg, 73.20.-r, 85.75.-d 

Recently, a theoretical work has played a key role in 
discovery of novel topological insulators (TI), where time 
reversal symmetry protects unusual surface states which 
could lead to electronic devices with new functionalities 
[TH3]. Unlike other fields of condensed matter physics, 
where most compounds are found haphazardly, here a 
number of new materials has been suggested based on 
band structure calculations [IHS], some of which have 
been confirmed experimentally [71 Hj . 

In addition to these unusual electronic properties, in- 
troducing magnetic order on the surface of TIs can lead 
to a gap for surface states, leading to completely insu- 
lating behavior. Here, novel phenomena such as the half 
quantum Hall effect [S] and quantized magnetoelectric 
effect [To] appear. The latter is parameterized by the 
quantized value of the Axion electrodynamics parame- 
ter 9 = TT, [101 lllj . When time reversal symmetry, as 
well as spatial symmetries like inversion are broken in 
the bulk, the parameter is no longer quantizcd[12, 13]. 
However, in most known materials, it is extremely small 
eg. = 10~^, 10""' in Cr203 and BiFeOs respectively 
[IB] . Thus, achieving large values of the parameter, 
even in the absence of quantization, is an important ma- 
terials challenge. For sufficiently large 0, > 7r/2, it has 
been proposed that magnetic domain walls will be associ- 
ated with conducting channels at the sample surface [l6]. 
Recent experimental progress in introducing magnetism 
via dopants [T31[Tn] in Bi2Se3 has been reported, as well 
as tuning of chemical potential into the small gap in- 
duced on surface states fTSI. Ideally however, magnetic 
topological insulators would combine band topology with 
intrinsic magnetic order, leading to large surface energy 
gaps. Unfortunately, most of known/proposed TIs in- 
volve p-electron orbitals, whose Coulomb interaction is 
weak and cannot support magnetism. Elements with Sd 
and Ad electrons do have large electronic correlations, but 
only small SOC. While both strong Coulomb interaction 
and large SOC is found in 4/ and 5/ electron systems, 
they usually form narrow energy bands, thus making it 



also hard to realize TI state of matter (see [T7] for an 
exception) . 

Due to the interplay of strong SOC and moderate elec- 
tronic correlations, very recently systems with 5d elec- 
trons have received a lot of research attention |18H23| . It 
has been suggested that pyrochlore iridates exhibit novel 
phases such, e.g., TI[2Tll22|, and topological semi-metal 
(TSM) '23] behavior. As in pyrochlores ^2-^207, spinel 
compounds AB2O4 have their four _B-sites forming a cor- 
ner sharing tetrahedral network. Furthermore, it may 
be expected that 5d systems in the spinel structure will 
be more tunable by pressure, external fields or by dop- 
ing as compared to the closely packed pyrochlore lattice. 
Previous studies suggested that electronic d^ configura- 
tion can realize an insulating state in 5d compounds [T^]. 
Based on valence arguments we therefore focus on spinel 
osmates ^Os204 {A is alkali metal element such as Mg, 
Ca, Sr, Ba), and investigate their electronic structure and 
magnetic properties using theoretical density functional 
calculations. Our numerical results reveal that the spinel 
structure can at least be a metastable crystallographic 
state for 5d elements which is additionally supported by 
a very recent synthesis of an empty iridium spinel com- 
pound Ir204[24]. Most interestingly, we find that the 
ground state is ferromagnetic and the electronic proper- 
ties are sensitive to the ^-site. Depending on Coulomb 
correlation, they show some exotic properties, and for the 
U relevant to osmate, we predict Ca0s2 04 and SrOs204 
behave as magnetic axion insulators with = tt. The 
quantized value is protected by the inversion symmetry 
of the lattice. However, since the surfaces break inver- 
sion, typically they are gapped allowing us to define 0. 
Unlike other proposals for generating large 0, which ex- 
ploit proximity to a topological insulator state [T5], here 
a generalized parity criterion [551 US] allows us to obtain 
large even when no proximate TI state exits. 

We perform our electronic structure calculations for 
the spinel Osmates based on local spin density approxi- 
mation (LSDA) to density functional theory (DFT) with 
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FIG. 1: Fragment of spinel crystal structure AOS2O4 (only 
Os and O atoms are shown) with the Os-O-Os angle being 
optimized. 



the full-potential, all-electron, linear-muffin-tin-orbital 
(LMTO) method[27,. Despite the 5d orbitals are spa- 
tially extended, recent theoretical and experimental work 
has given the evidence on the importance of Coulomb in- 
teractions in 5d compounds [1 9 . We utilize LSDA+U 
scheme [5S] to take into account the effect of Coulomb re- 
pulsion, and vary parameter U between and 1.5 eV. We 
use a 24x24x24 k-mesh to perform Brillouin zone inte- 
gration, and switch off symmetry operations in order to 
minimize possible numerical errors in studies of various 
(non-)collinear configurations. As the experimental lat- 
tice parameters are not available, we search for the stable 
crystal structures by locating the minimum in the calcu- 
lated total energy as a function of the lattice constant 
and internal atomic coordinates. 

Spinel structure (see Fig. 1) forms space group FdSm. 
In order to allow its relaxation, we change the lattice 
constant from 13 a.u. to 18 a.u. and vary the Os-O-Os 
bond angle from 90" to 120*^ with a step of l''. We con- 
firm that the choice of Coulomb U has only small effect 
on the determination of lattice parameters, and list the 
lowest energy structures obtained by LDA-f-SO (U=0) 
calculations in Table I. For comparison, we also list the 
same parameters for the pyrochlore iridate Y2lr207. We 
see that the ^-site element has a considerable effect on 
the Os-0 bond angle and its length. As we discuss be- 
low, this allows us to control the electronic structure in 
Osmates and design exotic topological phases. 

The 3s band of Mg is lower in energy, and appears 
around the Fermi level, which makes MgOs204 always 
metallic. We therefore do not discuss this compound here 
and concentrate our study on CaOs2 04. Its band struc- 
ture from non-magnetic LDA-I-SO calculation is found 
to be metallic and shown in Fig. 2a. The energy bands 
around the Fermi level appear as Je//=l/2 states similar 
to the ones found in Sr2lr04[Tn], and also in ¥911967 ^5]. 
where a metal rather than the interesting topological in- 
sulator scenario of [21: was obtained due to a 2-4-2 se- 



TABLE I: Theoretically determined structure parameters of 
spinel osmates. Angle denotes the Os-O-Os bond angle; Os- 
O and Os-Os denote the nearest neighbor Os-O bond length, 
and Os-Os bond length, respectively. 





Angle 


Os-Os Os-0 


CaOs2 04 


98.5" 


3.07 A 2.04 A 


SrOs2 04 


94.6° 


3.12 A 2mA 


BaOs204 


103.1" 


3.28 A 2.nA 


Y2lr204 


129.7" 


3.60 A I.99A 



quence of degeneracies at the T point. As importance of 
electronic correlations for 5d orbitals has been recently 
emphasized [H], we therefore perform LSDA+U+SO cal- 
culations. Although, the accurate value of U is not 
known, the Os-Os bond length of spinel osmates is 
shorter than that of Y2lr2 07, and one can expect that 
the U in CaOs2 04 is smaller than in Y2lr2 07. We there- 
fore believe that the U is in the range between 0.5 and 1.5 
eV. As in the pyrochlore structure, the Os spinel sublat- 
tice is geometrically frustrated. Naively one may expect 
that the magnetic configuration of CaOs2 04 is also non- 
conlinear as recently found in Y2lr207[23]. To search 
for possible magnetic ground states, we perform calcu- 
lations by starting with a number of different conlinear 
and noncolinear magnetic configurations including ferro- 
and antiferromagnetic (FM/AFM) colhnear (010), (110), 
(111), as well as non-collinear all-in/all-out, 2-in/2-out, 
3-in/l-out and some perpendicular configurations pro- 
moted by Dzyaloshinsky-Moriya interactions p51. We 
find that when U is less than 1.3 eV, only the FM-(OIO)- 
configuration retains its initial input magnetization di- 
rection; in all other configurations the moments depart 
from their input orientation. We also consider a two-up, 
two-down state which is suggested by the strong coupling 
limit where Os-0 bonds are nearly 90°. As shown in ^50], 
this leads to a 'Kitaev type' ferromagnetic interaction for 
spin components perpendicular to the plane. Although 
this structure is found to be stable for U < 0.8eV, it is 
higher in energy than the ferromagnetic state as shown 
in Table II. Regardless of the value of U, the FM con- 
figuration with magnetization along (010) is found to be 
the ground state, and the energy difference between this 
and other configurations is quite large. 

This is understood by examining lattice parameters in 
Table I, where the main difference between pyrochlore 
iridates and spinel osmates is the Os-O-Os bond an- 
gle and the Os-Os bond length. Due to the extended 
nature of bd orbitals, the 5d-2p hybridization is strong 
and important for the inter-atomic exchange interaction. 
In Y2lr207, the Ir-O-Ir bond angle is much larger than 
90" and Ir-O-Ir antiferromagnetic (AFM) superexchange 
interaction is dominant. This results in a strong mag- 
netic frustration and non-conlinear ground state mag- 
netic configuration [23j . In contrast to pyrochlores, the 
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TABLE II: The spin (S) and orbital (O) moment (in /ig) 
as well as the total energy Etot per unit cell (in meV) for 
several selected magnetic configurations of CaOs2 04 as cal- 
culated using LSDA+U+SO method with U=0.5 eV. {Etot is 
defined relative to the ground state.) The IDM is a co-planar 
configuration predicted for one sign of Dzyaloshinski Morya 
interactions in Ref. [29] 

Configuration: (010) all-in/out 2-in/2-out 2-up/2-dn IDM 



0.59 


0.30 


0.40 


0.25 


0.35 


0.14 


0.21 


0.22 


0.14 


0.16 


0.00 


290 


280 


518 


351 



TABLE III: Calculated parities of states at Time Reversal 
Invariant Momenta (TRIMs) of CaOs2 04. Only the 4 empty 
t2g bands are shown in order of increasing energy. Lx 2-K/a{- 
0.5,0.5,0.5), Ly 27r/a(0.5,-0.5,0.5), Lz 27r/a(-0.5,0.5,0.5) and 
L 27r/a(0.5,0.5,0.5)). 

Lx Ly Lz L 

U=0.5 eV ++++ ++++ +++H 

U=1.5 eV -+++ -+++ -+++ H 



pied states at the TRIM points i with odd parity. This 
allows us to analyze their topological electronic proper- 
ties based on the parities of occupied bands. One must 
also verify the Chern number of the bands vanish, to de- 
fine the magneto-electric coupling, which can be done by 
studying the evolution of the band structure from a triv- 
ial insulating state. The local environment of Os ion in 
spinels is oxygen octahedron, and the crystal-field split- 
ting of Os 5d states is large, which makes the e^ band 
higher (about 2 eV) than the Fermi level. The bands 
close to the Fermi level are mainly contributed by t2g 
mixed with O 2p orbitals. There are 4 Os per unit cell, 
thus the number of t2g bands is 24. If all of 24 t2g bands 
are fully occupied, the system cannot possess any topo- 
logically non-trivial properties. Noticing that Os occurs 
in its 3^ valence, for an insulator to occur, there are 20 
occupied and 4 empty t2g bands. Thus, instead of an- 
alyzing 20 occupied bands, we alternatively look at the 
parities of 4 empty bands. 

For the spinel lattice, the eight TRIM points are F, 
X, Y, Z, and four L points. For the F point, all of the 
t2g bands possess even parity, while for X, Y, and Z 
points, two of the 4 empty bands have even parity and 
the other two bands possess odd parity. Thus, regardless 
the value of (7, the summation of the number of empty t2g 
states with odd parities at F, X, Y, and Z points is 6. On 
the other hand, the Coulomb interaction has a significant 
effect on the parities of the bands around the L points. 
This is shown in Table III, where the summation of the 
number of the odd states are 6 and 4 for U=1.5 and 
0.5 eV, respectively. Therefore, according to Eq{T] for 
U=1.5 eV, we have — 0, which corresponds to a normal 
insulator. However, for U=0.5 eV, we have 9 — it, which 
leads us to the axion insulator with novel magnetoelectric 
properties. 

Since topological insulators must be separated from 
trivial insulators by a semi-metallic state [23], with 3D 
Dirac like dispersion, we additionally perform calcula- 
tions for a number of intermediate values of U to find out 
the boundary between semi-metal and Axion insulator as 
well as between semi-metal and Mott insulator. Our cal- 
culations for U's varying from 1.0 to 1.3 eV show that 
there are 3D Dirac crossings close to (0.02, x, 0.02)27r/a 
points of the BZ. where the value of x will change with 
U. With both decreasing and increasing the U, the Dirac 
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FIG. 2: Electronic band structure of CaOs204 shown along 
high symmetry direction, (a) LDA+SO; (b) LSDA+SO-f U, 
with U=0.5 eV; (c) LSDA+SO+U, with U=1.5 eV 



Os-O-Os angle in spinel osmates is close to 90*^, while 
Os 5d orbitals have stronger overlap, making a direct 
ferromagnetic exchange to be important for CaOs204. 

For the values of U=0.5 and 1.5 eV, electronic band 
structures along high symmetry lines appear to be insu- 
lating as shown in Fig. 2b and Fig. 2c, respectively. To 
check a possibility of the gap closure and the metallic 
behavior away from the high symmetry lines of BZ, we 
perform the calculation with 100 x 100 x 100 k-mesh. This 
very dense k-mesh confirms that both U=0.5 and 1.5 eV 
calculations show the band gaps of 0.01 and 0.08 eV re- 
spectively. 

In order to examine whether the system becomes an 
Axion insulator, we recall that similar to the Fu and 
Kane criterion [5], which is designed for non-magnetic 
systems with inversion symmetry, one can still use the 
parity eigenvalues for the states at eight Time Reversal 
Invariant Momenta (TRIM) in order to classify the mag- 
netic insulators with inversion symmetry and evaluate 
the magneto-electric coupling [251 [^ . 

9 = TT- Af (mod 2) (1) 

where Af — (X^fc^i)/^; Ni is the number of occu- 
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FIG. 3: Sketch of the predicted phase diagram for spinel os- 
mates. 



point moves and annihilates by meeting with other Dirac 
points, thus opening the energy gap and forming either 
the 9 ~ Ti Axion insulator or normal insulator, respec- 
tively. To summarize, our electronic phase diagram (see 
Fig. 3) for U <0.4 eV, predicts CaOs204 to be a metal; 
for 0.4< U <0.9 eV, an axion insulator; for 0.9< U <1.4 
eV, a topological semi-metal metal; for U >1.4 eV, a 
Mott insulator. Note, unlike in the cubic pyrochlore iri- 
dates [23], where the topological semi-metal has Dirac 
points exactly at the Fermi energy, here, because of lower 
symmetry of the magnetic order, the Dirac points are 
slightly above or below the Fermi energy. To see the 
sensitivity of these results to the lattice parameters we 
perform calculations with increasing and decreasing the 
volume by 6%, as well as adjusting the internal coordi- 
nates and performing small changes in Os-O-Os bond an- 
gle. We confirm that while there are some little changes 
in calculated energy bands, our predictions on Axion in- 
sulator behavior are robust. 

Replacing Ca by Sr, both the Os-O-Os bond angle and 
the bond length will change as shown in Table I. However, 
same as CaOs2 04, SrOs2 04 shows the same rich phase 
diagram as a function of U. 

We finally study BaOs204 which has both the largest 
bond angle and bond length as shown in Table I. These 
differences significantly affect its band structure: the 
LDA-f SO calculation with U=0 gives that at the F point 
the eight Jgyy = 1/2 states have degeneracies 4-2-2 and 
not 2-4-2 as found in Ca and Sr cases. Note that this se- 
quence of levels is the same as recently suggested by Pesin 
and Balents to realize topological insulating scenario in 
pyrochlore iridates [21]. Unfortunately, the bands cross- 
ing the Fermi level exist and cannot be removed by slight 
adjusting the lattice constant. Same with CaOs204 and 
SrOs2 04, considering the Coulomb interaction U will in- 
duce magnetism, but we do not find axion-insulator or 
Dirac- metal state for BaOs2 04 for any reasonable U. 

To summarize, using density functional based elec- 
tronic structure calculations, we have found that depend- 
ing on the strength of the Coulomb correlation among 
the 5d orbitals several exotic electronic phases may be 
realized in Os based spinel compounds. In particular a 



magnetic topological insulating phase (axion insulator) 
may be realized in CaOs204 and SrOs204 with large 
orbital magneto-electric parameter 9 — n. This re- 
search suggests that new functionalities such as control- 
ling electrical ^6J and optical [3Ij properties via magnetic 
textures, can be found in the new 5d spinel materials that 
we propose. 
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